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AN EXPLORATORY STUDY OF PARALLEL-STAGE SEPARATION 


OF REUSABLE LAUNCH VEHICLES* 

By John P. Decker and Joseph Gera 
Langley Research Center 

SUMMARY 

An exploratory study has been made of the stage separation of parallel- staged 
reusable launch vehicles. Static longitudinal aerodynamic data were obtained for both 
stages of a representative two-stage rocket-powered reusable-launch-vehicle concept 
when the stages were in close proximity to each other. The effects of vertical spacing, 
longitudinal spacing, and incidence angle were determined at Mach numbers of 3 and 6. 

In an attempt to identify the problems associated with mutual stage interference during 
the separation maneuver, the equations of longitudinal motion were numerically inte- 
grated. The calculation was carried out for a limited number of initial conditions by 
using a portion of the wind-tunnel data along with estimated values of the dynamic 
derivatives. 

The experimental results indicated large interference increments on the static - 
stability and normal-force characteristics of both vehicle stages. The static aerody- 
namic data for both vehicle stages were found to be dependent on the position and attitude 
of the second stage, the Mach number, and the relative sizes of the two stages. The 
relative sizes of the stages are dependent on the launch mode, the mission requirements, 
and the choice of propulsion system for the launch vehicle. 

The trajectory results indicated a strong dependence of the dynamic derivatives, 
initial attitude, and dynamic pressure on the separation maneuver. These results also 
indicated that a potentially hazardous situation could be expected when parallel-arranged 
stages are separated. However, to evaluate the potential of safely separating any given 
system will require a critical staging analysis of the system under consideration to 
determine the constraints for which safe separation may be achievable. Furthermore, 
other avenues of approach, such as trajectory shaping, vehicle shaping, various thrusting 
maneuvers including attitude control thrust, stability augmentation, and other auxiliary 
devices need to be evaluated. 

*The information presented herein was also presented in a thesis in partial fulfill- 
ment of the requirements for the degree of Master of Aerospace Engineering, University 
of Virginia, Charlottesville, Virginia, 1968. 


INTRODUCTION 


Interest in the potential of the recovery and reuse of launch vehicles during the 
past several years has resulted in a number of studies of various types of recoverable 
systems. The spectrum of concepts has varied from recoverable ballistic to winged 
reusable airbreathing vehicles. This paper is restricted to the class of reusable launch 
systems where the stages are arranged in parallel for launch. If staging must occur 
within the sensible atmosphere for these systems, aerodynamic interferences between 
the two vehicles may result in significant effects on the behavior of each vehicle during 
the staging maneuver. The Langley Research Center has therefore undertaken an 
exploratory investigation to ascertain these interferences and to interpret their influence 
on the staging maneuver. The discussion of reusable launch vehicles of reference 1 also 
indicated the existence of a potential problem in separating parallel -arranged stages of 
a reusable launch vehicle system. The present paper is an extension of the work pre- 
sented in reference 2 and is an initial attempt to explore the complex problem of 
parallel- stage separation. 

In order to provide meaningful information on the magnitude and character of the 
staging problem, it is necessary to obtain aerodynamic test data on each of the two vehi- 
cles in close proximity at conditions which might occur during staging and then to esti- 
mate the effects of these forces and moments on the relative motion of the two vehicles. 
For the present investigation, which is primarily exploratory, static longitudinal data 
were obtained on two representative stages of a rocket- powered reusable-launch-vehicle 
system when the vehicles were in close proximity to each other to obtain some under- 
standing of the aerodynamic-interference phenomena. The first stage was a simplified 
wing-body configuration whereas the second stage was a lifting body. Vertical and longi- 
tudinal spacing, as well as incidence angle, were varied for a range of angle of attack and 
at Mach numbers of 3 and 6. These data were then used in integrating the equations of 
motion of both vehicles to determine the relative behavior of both stages during a staging 
maneuver. This analysis was made by using the coupled longitudinal equations of motion 
for both stages. In these equations, the general load terms were replaced by the experi- 
mental aerodynamic data and estimates of the damping derivatives. For selected values 
of stage characteristics and initial conditions, the equations were solved numerically. 

SYMBOLS 

The aerodynamic loads of each stage have been resolved and reduced into coeffi- 
cient form with respect to its own body-axis system. The location of the origin, and the 
orientation of the body axes are shown for the two body-axis systems in figure 1. The 
geometric constants used in reducing the aerodynamic loads into coefficient form are 
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given in table I. The physical quantities defined in this paper are given in both the U.S. 
Customary System of Units and the International System of Units (SI). Conversion fac- 
tors between the two systems are given in reference 3. 
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planform area of second- stage pitch control flaps, feet2 (meter s2) 

axial-force coefficient, Measured axial force 

| pM 2 S 


pitching- moment coefficient, Measured pitching moment 

| pM 2 Sc 

damping in pitch, 
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rate of change of pitching moment due to acceleration along Z-axis, 
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in, seconds2/foot (seconds2/meter) 


normal-force coefficient, 


Measured normal force 
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rotary stability derivative, — — 

8 £L£L 

2V 

normal-force slope, per degree 


thrust coefficient, 


Thrust force 
^pM 2 S 


reference length (see table I); for the first stage it is equal to the mean aero- 
dynamic chord, for the second stage it denotes the body length, feet 
(meters) 

distance between centers of mass, feet (meters) 

perpendicular distance from center of mass to line of thrust, feet (meters) 


F X ,F Z components of resultant force, pounds (newtons) 
g acceleration due to gravity, feet/second2 (meter s/second2) 


h altitude, -z e , feet (meters) 
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mass moment of inertia, slug-feet 2 (kilogram -meters 2 ) 


i incidence angle, " a l> degrees 

k ratio of specific heats 

l body length, feet (meters) 

M Mach number 

My total pitching moment, foot-pounds (meter-newtons) 

m mass, slugs (kilograms) 

p static pressure, pounds/foot2 (newtons/meter 2 ) 

q angular velocity, radians/second 

R radius, inches (centimeters) 

S reference area (see table I), feet 2 (meters 2 ) 

t time, seconds 

u,w components of velocity along body axes, feet/second (meters/second) 

V resultant velocity, (u 2 + w 2 )^^, feet/second (meters/second) 

X,Y,Z coordinate axes 

x,z distances along body axes, feet (meters) 

x e ,z e distances along earth-fixed axes, feet (meters) 

x a ,xjj,r a ,z a fuselage coordinates (see fig. 1(b)), inches (centimeters) 

Ax,Az separation variables (positive when center of gravity of second stage is 

ahead and above center of gravity of first stage, see fig. 4), feet (meters) 

a angle of attack, degrees 
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Y flight-path angle, 9 - a, degrees 

6 e resultant angle of pitch control flap (positive when trailing edge is down), 

6e, right + 6e,left > degrees 

e angle between X-axis and thrust line, degrees 

9 angle of pitch, degrees 

p atmospheric density, slugs/foot3 (kilograms/ meter3) 

At time increment, seconds 

A dot over a symbol denotes first derivative with respect to time. 

The subscripts 1 and 2 indicate whether the physical quantity represented by the 
principal symbol refers to stage 1 or stage 2. 

APPARATUS AND TESTS 

Models 

Details of the launch- vehicle configuration are shown in figure 1. The launch 
vehicle consisted of a simplified wing-body first stage with a lifting-body second stage. 
The second stage was placed on top of the first stage, the flat bottom being parallel with 
the first- stage wing upper surface. The longitudinal location of the second stage was 
such that its moment reference center was behind the first- stage moment reference cen- 
ter in the normal carrying position, that is, A x/l} = -0.051. Photographs of the two- 
stage system prior to staging and at staging conditions are shown in figure 2. 

The first stage consisted of a semicylindrical fuselage with an ogival forebody and 
delta wing. The wing had 65° leading-edge sweep and the airfoil section was a half- 
diamond section that had a maximum thickness of 4 percent of the local chord at the 
40-percent-chord station. The wing was flat on the upper surface, and thus the wing had 
negative camber. No longitudinal or vertical control surfaces were provided on the first 
stage for the present tests. 

The second stage was a lifting-body type of vehicle with cross-sectional shape pro- 
gressions along the X-axis as shown in figure 1(c). The second stage had essentially a 
flat bottom and was provided with vertical stabilizing surfaces. Pitch controls were pro- 
vided on the second stage which consisted of upper and lower surface flaps located near 
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the base of the body. (See fig. 1(c).) The ratio of planform area of the pitch controls to 
the total planform area was 0.113 for the upper surface controls and 0.187 for the lower 
surface controls. 


Support Mechanism 

Separate sting supports were provided for the first and second stage, preset ver- 
tical movement between the stages being provided in the vertical plane by the support 
system to which the stings were attached. (See fig. 3.) Longitudinal movement and 
incidence angle between the stages were provided by using spacers and sting adapters on 
the upper sting support. The complete support apparatus was attached to an arc strut 
which varied the angle of attack of the vehicles. 

Tests 

The wind-tunnel tests were conducted at nominal Mach numbers of 3 and 6 in the 
2-foot hypersonic facility at the Langley Research Center described in reference 4. 
Position variables Ax and Az, as well as the relative incidence angle i, were varied 
for an angle- of- attack range of approximately -9° to 12°. (See fig. 4.) Static longitu- 
dinal aerodynamic force and moment data were simultaneously obtained for the first 
and second stage by use of individual internal six-component strain-gage balances. No 
composite configurations, that is, with the first stage and second stage connected, were 
tested. 

All data were obtained with the model smooth; that is, no boundary-layer transition 
strips were used. At the Reynolds numbers of these tests, laminar flow may be expected 
to exist over most of the model. Individual vehicle angles of attack were corrected for 
balance and sting deflection under load. No base drag corrections were made for either 
the first or second stage. The average tests conditions and Reynolds number variations 
are as follows: 


Mach 

Stagnation pressure 

Stagnation temperature 

Reynolds number 
based on a 1-foot 
(0.3048 m) length 

number 

atm 

kN/m^ 

°R 

°K 

3 

0.5 

50.6 

560 

322 

0.8 x 10 6 

6 

3.0 

303.6 

760 

422 

1.0 


METHOD OF ANALYSIS 

The separation maneuver of two parallel stages is illustrated schematically in fig- 
ure 5. At release, the second stage would be at an initial spacing distance and attitude 
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with respect to the first stage. A trapeze or similar mechanism could be employed to 
achieve the desired release conditions. Potential problems that may arise during the 
staging sequence are illustrated by the lower two sketches of figure 5. The divergence 
of the center of gravity alone would not imply safe separation since the second stage may 
rotate into the first stage. Consequently, realistic analysis can only be accomplished if 
the separation maneuver is regarded as the motion of rigid bodies. 

Two parallel lifting stages, separating from each other aerodynamically, represent 
a complex dynamic system. For the complete description of this system, not less than 
12 degrees of freedom must be taken into account. Although this is possible, at least 
theoretically, the final accuracy of the analysis will depend on the correct analytical 
representation of the aerodynamic forces and moments. 

From the available data, the aerodynamic loads cannot be predicted if any lateral 
motion occurs during the separation maneuver; therefore, it is assumed that the motion 
of both stages takes place in the common plane of symmetry. With this restriction, the 
Euler equations of motion (ref. 5) with respect to the body axes for either stage are 

m(u + wq) = F x - mg sin 9 

m(w - uq) = Fg + mg cos 6 \ (1) 

Iyq= My J 

The angular orientation and position of either stage are expressed as 

e = q "1 

x e = u cos 9 + w sin 9 \ (2) 

z e = -u sin 9 + w cos 9J 

This set of six equations is completely general and valid even for large distur- 
bances. However, in the development of the expressions for the force and moment 
terms, the following assumptions are made: 

(1) Since the time lapse between the initiation of the separation maneuver and full 
separation (or collision) of the two stages is small, the response of the system to any 
control deflection is assumed to be negligible. Hence, the controls are considered to be 
fixed. 

(2) During separation, the effects of speed changes on the aerodynamic-force 
coefficients are neglected. 

(3) Cross coupling between the two stages takes place only through the normal- 
force and pitching- moment coefficients. Consequently, the axial-force coefficient of 
either stage was assumed not to be affected by the proximity of the other stage. 
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With these assumptions, the force and moment terms may be expressed as 
F X = |pV 2 s(c T cos e - C A ) 

F Z = ipV 2 s|c T sin e - C N + JL C Nq q ) > (3) 

My = |pV 2 Sc(^ C T + C m + ^ C mq q + C m ^ 

In equations (3) the rotary stability derivatives Cjq q and C mq and the derivative 
Cm^ are considered to be constants for each stage during the separation maneuver since 
these derivatives are not presently available for two vehicles in close proximity. How- 
ever, the influence of these derivatives on the separation maneuver can be investigated 
by varying the magnitude of these quantities. Orientation of the thrust vector is specified 
by e (the angle between the X-axis and the thrust vector) and dp (the perpendicular 
distance from the center of gravity to the thrust line). These quantities, as well as the 
thrust coefficient C^, were assumed to be constant for each stage since the time lapse 
between the initiation of the separation maneuver and full separation (or collision) of the 
two stages is small. Equations (3) are arbitrary to a large extent, and they were written 
with the experimental part of this investigation in mind. The aerodynamic coefficients, 

C A , Cjvj, and C m , represent the experimental data measured on each stage. 

Up to this point, the form of the governing equations has been identical for both 
stages. For each stage, however, the aerodynamic coefficients are defined separately, 
that is, 

C A,1 = C A , !(<*!) 

C A,2 = c A,2(“l) 

C N,1 = c N> l(“ 1 ,« 2 » A x,Az) 

, \) (4) 

C N,2 = c N,2( a l» Q! 2» Ax » Az ; 

C m ,l = C m> i(a 1 ,a2,Ax,Az) 

c m,2 = Cm, 2 ( a l*“2» Ax » Az ) 

From the experimental data of the present investigation, these relations can be 
constructed explicitly in the form of tabular functions. The separation variables, a p 
a 2 , Ax, and Az, are illustrated in figure 4. 

The resulting equations of motion were integrated numerically on a digital com- 
puter. The values of the aerodynamic coefficients were obtained by linear interpolation 
between the discrete points of the tabular functions at each integration step. The 
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atmospheric density p in equations (3) corresponds to that given by the U.S. Standard 
Atmosphere, 1962 (ref. 6) at each height. 

For the numerical integration, the size of the computing interval was carefully 
selected by the use of two different integration techniques for each of which the size of 
the computing interval was systematically varied. Comparison of the solutions obtained 
in this manner indicated that a value of At between the values of 0.005 and 0.01 second 
resulted in an acceptable level of accuracy. Consequently, a fourth-order fixed- interval 
integration technique was used. This integration technique uses the Adams-Bashforth 
predictor formula in conjunction with the Adams- Moulton correction formula. A fourth 
order Runge-Kutta formula with minimum error bounds was used to start the solution. 
(See refs. 7 and 8.) The integration did not proceed beyond t = 5 seconds since safe 
separation or collision occurred within this time. 

When the numerical solutions of the present analysis are interpreted, it should be 
remembered that in addition to the assumptions listed, equations (3) contain another 
source of potential error. Ideally, the unaccelerated flight of a single body can be com- 
pletely simulated in the wind tunnel, whereas the separating flight of two parallel stages 
is impossible to simulate either kinematically or dynamically with fixed models in the 
test section. Thus, the direct use of the measured coefficients in equations (3) involves 
the error due to the wind-tunnel simulation technique where the coefficient data could be 
obtained only at y^ = y 2 and Vj = V^- This error is lessened by using the magnitudes 
of the instantaneous velocities in equations (3). However, this procedure will not take 
into account the effects of different flight -path angles on the aerodynamic interference 
nor any scaling or Reynolds number effects. 

PRESENTATION OF RESULTS 

The results have been divided into two principal parts. These two parts consist of 
the experimental results obtained for the two stages in close proximity and the calculated 
results obtained by integrating the equations of motion for each of the two stages during 
a staging maneuver. Since the aerodynamic interferences that occur when two bodies 
are in close proximity are potential stability and control problems for both vehicle 
stages, the experimental results are used to discuss the implication of these interfer- 
ences on the stability characteristics of the stages when in close proximity, whereas the 
calculated results are used to discuss the potential effects on the vehicle behavior of the 
aerodynamic interference. 
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EXPERIMENTAL RESULTS 


The basic longitudinal aerodynamic characteristics of the first and second stages 
at different vertical and longitudinal spacings, incidence angles, and Mach numbers are 
presented in figures 6 to 9. For reference, the interference -free data, or the data at a 
very large separation distance (Az \/l-± = °°), are also shown in these figures. 

First Stage 

Second-stage m oment reference center aft of firsL-stage moment reference center; 
Ax/Zj = -0.051.- Figure 6(a) presents the effect of vertical spacing at M = 3, i = 0°, 

second-stage moment reference center aft of first-stage moment reference center, and 
= -0.051. This figure indicates that the static stability level did not change appre- 
ciably with vertical spacing. However, there is a large positive increment in C m \ due 
to the interference and a significant change in the angle for trim. Also, for the maximum 
vertical spacing of the tests, large interferences are still present because the curve at 
Az Jl\ = 0.224 has not approached the interference -free curve (Az/Zj = °o^. Figure 6(a) 
also shows that there are large negative increments in Cj^i due to the interference 
and no appreciable change in Increasing i from 0° to 5° (fig. 6(b)) and 10° 

(fig. 6(c)) at M = 3 resulted in further positive increments in C m i and further nega- 
tive increments in Cjj } i. (Compare C m i and Cjj j at Az /l± = 0.127, for example.^ 

At M = 6 and i = 0°, figure 7(a) shows that both the magnitude of C m j and 
stability level of the first stage varied with vertical spacing. This result is in contrast 
to the data at M = 3 (fig. 6(a)) where only the magnitude of C m> i was shown to vary 
with vertical spacing. Also in contrast to the data at M = 3, the data at the largest 
vertical spacing (az jl-^ = 0.224^ at M = 6 (fig. 7(a)) indicates that the first stage is 
approaching interference -free conditions since both the normal-force and pitching- 
moment curves are approaching the interference -free curves. However, at this vertical 
spacing and also at Az/Zj = 0.199, the first stage has pitch-up tendencies at a j greater 
than 3°. Increasing i from 0° to 10° (figs. 7(a), 7(b), and 7(c)) at M = 6 caused posi- 
tive increments in C m i and a degradation of the static stability at all vertical spacings. 
Figures 7(a), 7(b), and 7(c) also show that the presence of the second stage caused nega- 
tive increments in and some localized changes in 

interference-free conditions. 

Second-stage mo ment reference center forward of first-stage moment reference 
center; Ax/Zj = 0.160.- At the forward longitudinal spacing Ax/Zj = 0.160 and at 
M = 3 and i = 0° (fig. 6(d)) not only did the magnitude of C m j vary with vertical 
spacing, but also the stability level decreased with increasing vertical spacing so that at 
Az^Zj = 0.199 and 0.244, the first stage is statically unstable at less than about 2°. 


(<*«), as compared with 
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Figure 6(d) also shows that in comparison with the interference -free curve, both (CNq,)^ 
and the magnitude of C^l at a-y greater than about -8° decreased because of the 
presence of the second stage. However, increasing the vertical spacing did not appre- 
ciably change either (CNq,)^ or C^ Increasing i from 0° to 10° (figs. 6(d), 6(e), 

and 6(f)) caused further decreases in the stability level, the magnitude of C m and the 
magnitude of Cj^j, as compared with the interference -free curves. 

At M = 6 (figs. 7(d), 7(e), and 7(f)), the first stage was statically unstable at the 
smaller vertical spacings for less than about 2°. These same figures also show 
that the magnitude of C m ^, C^j ^, and varied with vertical spacing. 

Second Stage 

Second-stage moment reference center aft of first-stage moment re ference center; 
A x/Zj = -0.051.- For the second stage at M = 3, i = 0°, second-stage moment reference 
center aft of first-stage moment reference center, and A x/l\ = 0.051 (fig. 8(a)), both 
the stability level and magnitude of C m> i varied with vertical spacing. At Az/Zj 
less than 0.175, the second stage was approximately neutrally stable or statically unstable 
at all a\. Coupled with these changes in stability and magnitude of C m 2 are changes 
in both Cj^ 2 and (Cno^’ Increasing i from 0° to 10° at M = 3 (figs. 8(a), 8(b), 

and 8(c)) caused the second-stage stability level to decrease and caused large negative 
increments in the magnitude of C m> 2 and large positive increments in Com- 

parison of the normal-force and pitching-moment data (figs. 8(a), 8(b), and 8(c)) with the 
corresponding data for the first stage (figs. 6(a), 6(b), and 6(c)) indicates that increasing 
i would tend to cause the first stage to nose up but away (-Cn ? i) from the second stage, 
whereas the second stage would tend to nose down but away (+ 0 ^ 2 ) from the first stage. 

Comparison of the data for the second stage at M = 3 (figs. 8(a), 8(b), and 8(c)) 
with the data at M = 6 (figs. 9(a), 9(b), and 9(c)) at similar geometric conditions indi- 
cates that the results at both Mach numbers are very similar. As would be expected, 
increasing the Mach number from 3 to 6 decreased the actual magnitude of the force and 
moment coefficients. 

Second-stage moment reference center forward of first-stag e moment reference 
center; Ax/Zj = 0.160.- At M = 3, second-stage moment reference center forward of 

first-stage moment reference center and A y^/l\ = 0.160, figures 8(d), 8(e), and 8(f) indi- 
cate that the second stage was statically stable at almost all vertical spacings, incidence 
angles, and aj. However, the presence of the first stage caused increments in both 
Cm, 2 and Cjj 2> the increments becoming larger with increasing i. Again, similar 
results are shown at M = 6. (See figs. 9(d), 9(e), and 9(f).) 
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Second- stage contro l effectiveness.- Because of the large interference increments 
found for the second stage, an investigation was made to determine the control effective- 
ness of the second stage at interference-free conditions. Figure 10 presents this control 
effectiveness data for control deflections of 0° and ±25°. Comparison of these data with 
the data in figures 8 and 9 indicates that there appears to be no reasonable control 
deflection which could overcome all the large pitching- moment increments shown for the 
second stage when in proximity to the first stage. 

Comments on Interference Increments 

Further analysis of the data shown in figures 6 to 9, along with schlieren data, has 
led to two major conclusions as to the source of the behavior of the force and moment 
coefficients as compared with the interference-free values. These conclusions are illus- 
trated in figure 11. The first conclusion is that the effects on the first stage are caused 
by the impingement of the principal disturbance generated by the second stage. Conse- 
quently, the changes in forces and moments previously shown for the first stage are 
approximately proportional to the strength of this disturbance and the area affected. The 
second conclusion is that the effects on the second stage are caused by the flow field from 
the first stage to which must be added the effects of the first reflection of the principal 
disturbance caused by the second stage. Furthermore, not only are the static aerody- 
namic data dependent on the relative position and attitude of the second stage and Mach 
number, but the data are also dependent on the relative sizes of the two stages since the 
flow field of the first stage and the strength of the disturbance generated by the second 
stage are a direct function of the sizes of the two vehicles. The relative sizes of the two 
stages are dependent on the launch mode, the mission requirement, and the choice of 
propulsion system for the launch vehicle. 

TRAJECTORY RESULTS 

The time involved in obtaining solutions of the equations of motion (eqs. (1) and (2)) 
for each stage put practical limits on the number of variables and vehicle characteristics 
investigated. Accordingly, staging was assumed to occur in climbing flight and the Mach 
number selected for the analysis was 3 since this value could correspond to either an 
abort or mission- staging condition. These considerations led to the selection of a num- 
ber of fixed vehicle characteristics and initial conditions given in table II. 

Presented in figure 12 is a typical altitude Mach number ascent trajectory plot 
shown only out to a Mach number of 8 and indicates where abort and mission- staging con- 
ditions might be required. Reusable launch vehicles using airbreathing engines may 
stage at the lower altitudes and consequently higher dynamic pressures whereas rocket- 
powered vehicles would stage at the higher altitudes and lower dynamic pressures. 
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Reference 1 has indicated that airbreathing vehicles may be able to stage at the lower 
dynamic pressures with a negligible performance decrement by executing a pull-up 
maneuver. Normal mission staging could occur, depending on the vehicle concept, any- 
where above a Mach number of about 3, whereas abort staging could occur at any Mach 
number below mission staging and at any altitude. 

Example Solutions 

Figure 13 illustrates typical numerical results at an altitude of 70 000 feet 
(21 336 meters) for two particular choices of the initial conditions and the damping-in- 
pitch characteristics. The vehicles were assumed to have safely separated when the 
vehicles were one first- stage body length apart and were assumed to have collided when 
the extremities of the vehicles touched each other. The main points of this figure are 
that (1) the time it took the vehicles to collide or separate was of the order of 3 seconds, 
a very short time for aerodynamic controls to be effective; (2) the net change in the ini- 
tial and final velocities of either stage during the separation maneuver would amount to 
a net Mach number change of about 0.10 and would have a negligible effect on the static 
aerodynamic coefficients; and (3) the difference between yj and y 2 which was a maxi- 
mum of 1.75° for the safe separation case, could introduce some errors in the calcula- 
tions since the static aerodynamic coefficients were obtained for yj = y<^ and could 
affect whether the vehicles would collide or safely separate. 

As was previously implied, no attempt was made in the present study to investi- 
gate the effect of all important vehicle characteristics, initial conditions, and damping 
characteristics either separately or in all reasonable combinations. Instead, attention 
was given to some of the least accurately estimated damping variables, namely, ^C m q^ 

and and to some of the initial conditions which have obvious bearing on suc- 

cessful separation. Included in the latter category were the initial relative attitude, 
dynamic pressure, and vertical thrusting. 


Effect of Dynamic Derivatives 


Figure 14 shows the effect of the dynamic derivatives where (c m q^ is plotted 
against (c m q^- Indicated in this figure are the regions where the two vehicles collided, 
an approximate region where they safely separated, and a region of uncertainty as to 
whether the two vehicles would collide or separate. Also indicated in this figure are the 
best current estimates of the interference-free values of C m q for both the first and 
second stages. These values are one or more orders of magnitude smaller than the val- 
ues needed to achieve safe separation at these conditions. Although the magnitude of 
Cn% and Cftjq for both the first and second stages have not been perturbed in this anal- 
ysis, it can be expected that similar figures would be generated for these quantities. 
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Effect of Initial Attitude 


For the purpose of illustrating some of the other important variables that may need 
to be considered to achieve safe separation, figure 15 shows the type of results that can 
be generated by selecting a value of C m q which lies in the safe separation region for 
both the first and second stages (square symbol, fig. 14). Here the incidence angle and 
the angle of attack of the first stage have been varied. Again, the regions of collision 
and safe separation are shown. 

Effect of Dynamic Pressure 

Figure 16 shows the effect of varying the dynamic pressure or altitude for the 
separation maneuver by taking a suitable safe separation value of i and a\ (circular 
symbol, fig. 15). The distance between the centers of gravity of the two stages in 2 sec- 
onds after release has been plotted against dynamic pressure. The dynamic -pressure 
range is for altitudes above and below the nominal altitude for a rocket-powered vehicle 
at M = 3. Two curves are presented; one is for the condition in which no thrust was 
used in the calculations and the other is for the condition in which a hypothetical vertical 
downward thrust corresponding to a thrust-weight ratio of 0.1 has been applied to the 
first stage. Using either vertical or longitudinal thrust would present jet-plume body 
interactions which, in these calculations, have not been evaluated. For the no-thrust 
condition and at zero dynamic pressure, the position and attitude of the vehicles do not 
change; and separation does not occur. Safe separation can be expected only at the 
higher dynamic pressures for the no-thrust condition. For the vertical-thrust condition, 
safe separation could result at low dynamic pressures but as the dynamic pressure 
increases, there is a region of collision. At the higher dynamic pressures, safe sepa- 
ration would also be predicted. 

Figures 14, 15, and 16 have indicated some of the problems associated with 
parallel- stage separation. These results have indicated that evaluating the potential of 
any parallel- staged system requires a critical dynamic analysis of the two vehicles when 
in close proximity to determine the constraints under which safe separation of the two 
vehicles may be achievable. 


CONCLUDING REMARKS 

An exploratory study has been made of the stage separation of parallel- staged 
reusable launch vehicles. Static longitudinal aerodynamic data were obtained for both 
stages of a representative two-stage rocket-powered reusable-launch-vehicle concept 
when the stages were in close proximity to each other. The effects of vertical spacing, 
longitudinal spacing, and incidence angle were determined at Mach numbers of 3 and 6. 


14 



In an attempt to identify the problems associated with mutual stage interference during 
the separation maneuver, the equations of longitudinal motion were numerically inte- 
grated. The calculation was carried out for a limited number of initial conditions by 
using a portion of the wind-tunnel data along with estimated values of the dynamic 
derivatives. 

The experimental results indicated large interference increments on the static- 
stability and normal-force characteristics of both vehicle stages. The static aerody- 
namic data for both vehicle stages were found to be dependent on the position and attitude 
of the second stage, the Mach number, and the relative sizes of the two stages. The rela- 
tive sizes of the stages are dependent on the launch mode, the mission requirements, and 
the choice of propulsion system for the launch vehicle. 

The trajectory results indicated a strong dependence of the dynamic derivatives, 
initial attitude, and dynamic pressure on the separation maneuver. These results also 
indicated that a potentially hazardous situation could be expected in separating parallel- 
arranged stages. However, to evaluate the potential of safely separating any given sys- 
tem will require a critical staging analysis of the system under consideration to deter- 
mine the constraints for which safe separation may be achievable. Furthermore, other 
avenues of approach, such as trajectory shaping, vehicle shaping, various thrusting 
maneuvers including attitude control thrust, stability augmentation, and other auxiliary 
devices need to be evaluated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 17, 1968, 

124-07-05-02-23. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF MODELS 
AND REFERENCE DIMENSIONS 


First stage: 

Body length, in. (cm) 23.760 (60.350) 

Overall length, in. (cm) 23.760 (60.350) 

Aspect ratio 1.865 

Span, in. (cm) 13.802 (35.057) 

Total wing area, in2 (cm2) 102.135 (658,934) 

Reference area, in2 (cm2) 102.135 (658.934) 

Root chord, in. (cm) 14.800 (37.592) 

Tip chord 0 

Mean aerodynamic chord, in. (cm) 9.867 (25.062) 

Reference length, in. (cm) 9.867 (25.062) 

Second stage: 

Body length, in. (cm) 11.465 (29.121) 

Overall length, in. (cm) 13.706 (34.813) 

Aspect ratio 0.678 

Span, in. (cm) 6.050 (15.367) 

Planform area less vertical fins, in2 (cm2) 53.980 (348.257) 

Reference area, in2 (cm2) 53.980 (348.257) 

Reference length, in. (cm) 11.465 (29.121) 
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TABLE n.- VEHICLE CHARACTERISTICS AND INITIAL CONDITIONS 


First stage: 

Weight, lb (N) 

Iy, slugs-ft2 (kg- m2) 

c, ft (m) 

S, ft2 (m2) 

C m ^, sec 2/ft (sec2/m) 

Cm q 

C Nq 

V^, ft/sec (m/sec) . . 

yi, deg 

l\, ft (m) 

Second stage: 

Weight, lb (N) .... 

Iy, slugs-ft2 (kg- m2) 

c, ft (m) 

S, ft2 (m2) 

C m ^, sec 2 /ft (sec2/m) 

C m q 

C N q 

Ax/lj 

Az A 

y 2 » de e 

Vj, ft/sec (m/sec) . . 


. . 345 000 (1 534 636) 

11.00x106 (14.9X106) 

78.11 (23.81) 

6401 (594.7) 

. . -0.00002 (-0.00006) 

Variable 

0 

3100 (945) 

26 

188.1 (57.33) 

, . 300 000 (1 347 811) 

1,75 xlO 6 (2.4x106) 

90.76 (27.66) 

3383 (314.3) 

. -0.000005 (-0.00001) 

Variable 

0 

-0.051 

0.127 

26 

3100 (945) 
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(c) Second stage. 
Figure 1.- Concluded. 
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(a) Prior to stage separation. 


<b) During stage separation maneuver. 

Figure 2.- Photographs of launch vehicle system. 
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variables. 
























(c) i = 10°: Ax/li = -0.051. 
Figure 7.- Continued. 
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(f) i = 10°; A x/lj = 0.160. 
Figure 7.- Concluded. 
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Figure 13.- Typical results 
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